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^ Abstract 

^ In this paper, using the Darboux transformation, we demonstrate the gener- 
ation of first order breather and higher order rogue waves from a generahzed 

I "I nonhnear Schrodinger equation with several higher order nonhnear effects 

C/^ representing femtosecond pulse propagation through nonlinear silica fibre. 

P5 The same nonlinear evolution equation can also describes the soliton type 

r-^ nonlinear excitations in classical Heisenberg spin chain. Such solutions have 

i_i a parameter 71 denoting the strength of the higher order effects. From the nu- 

_j merical plots of the rational solutions, the compression effects of the breather 

>■ and rogue waves produced by 71 are discussed in detail. 
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^ 1. Introduction 

> 

It is well known that one of the most challenging aspects of modern 

science and technology is the nonlinear nature of the system, which is con- 
sidered to be fundamental to the understanding of many natural phenomena. 
In recent years, nonlinear science has emerged as a powerful subject for ex- 
plaining the mysteries of the challenging nature. Nonlinearity is a fascinating 
occurrence of nature whose importance has been well appreciated for many 
years, in the context of large amplitude waves or high intensity laser pulses 
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observed in various fields ranging from fiuids to solid state, chemical, bio- 
logical, nonlinear optical and geological systems [IlEllSllllElElITlElin]. 
This fascinating subject has branched out in almost all areas of science, and 
its applications are percolating through the whole of science. In general, 
nonlinear phenomena are often modeled by nonlinear evolution equations 
exhibiting a wide range of high complexities in terms of different linear and 
nonlinear effects. In recent years, the advent of high-speed computers, many 
advanced mathematical software and development of many sophisticated and 
systematic analytical methods in the study of the nonlinear phenomena and 
also supported by many experiments have encouraged both theoretical and 
experimental research. In the last few decades or so, nonlinear science has 
experienced an explosive growth by the invention of several exciting and fas- 
cinating new concepts like solitons, dispersion managed solitons, dromions, 
rogue waves, similaritons, supercontinuum generation, complete integrabil- 
ity, fractals, chaos etc. [IlSISlSllSlElEIHlElIiniin]. Many of the completely 
integrable nonlinear partial differential equations (NPDEs) admit one of the 
most striking aspects of nonlinear phenomena called as soliton, which de- 
scribe soliton as a universal character and they are of great mathematical 
as well as physical interest too. The study of the solitons and other related 
issues of the construction of the solutions to a wide class of NPDEs have 
become one of the most exciting and extremely active areas of research in 
science and technology for many years. 

In addition to several developments in soliton theory, recent develop- 
ments in modulational instability (MI) have also been widely used to ex- 
plain why experiments involving white coherent light supercontinuum gener- 
ation(SCG), admit a triangular spectrum. Such universal triangular spectra 
can be well-described by the analytical expressions for the spectra of Akhme- 
diev breather solutions at the point of extreme compression. In the context 
of the NLS equation. Peregrine already in [12] had identified the role of MI in 
the formation of patterns resembling freak waves or rogue wave (RW); these 
theoretical results were later supported by several experiments. Rogue waves 
in ocean are localized large amplitude waves on a rough background, which 
have two remarkable characteristics: a) "appear from nowhere and disappear 
without a trace" ^13j, b) exhibit one dominant peak. RWs have recently ap- 
peared in several areas of science. Particularly in photonic crystal fibres, 
RWs have been well established in connection with SCG [14J. This actually 
has stimulated research for RWs in other physical systems and has paved the 
way for many important applications, including the control of RWs by means 



of SCG [151 [in], as well as studies in superfluid Helium [T^, Bose Einstein 
condensates [H] , plasmas [191 EO] , microwave [21] , capillary phenomena [22] , 
in telecommunication data streams [23], inhomogeneous media [21], water 
experiments [251 126] , and so on. More recently, Kibler et al. [27] using their 
elegant experimental apparatus in optical fibres were able to generate fem- 
tosecond pulses with strong temporal and spatial localization and near-ideal 
temporal Peregrine soliton characteristics. 

In the recent past, several equations have been shown to admit the rogue 
wave solutions. For example, in addition to the NLS equation, the Hirota 
equation [2H1 I2SI EO], the first type derivative NLS equation [21], the third 
type DNLS equation [32], the Fokas-Leneels equation [33], the NLS-MB equa- 
tions [31], the Hirota Maxwell-Bloch(MB) equations [35], the Sasa-Satsuma 
equation [36], the discrete Ablowitz-Ladik and Hirota equation [37], the two- 
component NLS equations [33 [321 HO], the three-components NLS equations 
[H], the variable coefficient NLS [1211131 IH iSl [16] , the variable coefficient 
derivative NLS [17], and the variable coefficient higher order NLS(VCHNLS) 
[IE], the rogue waves in dissipative systems |l9]are a few of the nonlinear evo- 
lution equations which admit RWs. From the above studies, it is clear that 
one of the possible generating mechanisms [50] for the higher order RW is 
the interaction of the multiple breathers possessing the same and the very 
particular frequency of the underlying equation. 

In recent years, there has been a considerable interest in the study of 
the nonlinear excitations of the spin chains with competing bilinear and bi- 
quadratic interactions. In particular, the complete integrability and nonlin- 
ear excitations of spin chains with spin magnitude S >1 has been established 
if suitable polynomials in {Si, Sj) are added to the original bilinear Heisen- 
berg spin Hamiltonian. In this connection and also from the mathematical 
point of view, it is of interest to study the influence of the biquadratic interac- 
tions on the integrability of the Heisenberg bilinear spin chain in the classical 
limit as well. Considering the above points, one of the authors of this pa- 
per, has investigated the integrability aspects of a classical one-dimensional 
isotropic biquadratic Heisenberg spin chain in its continuum limit up to or- 
der 0{a^) in the lattice parameter through a classical differential geometric 
approach [511 [521 ES] and investigated the soliton and integrability aspects 
of the corresponding generalized nonlinear Schrodinger Equation (GNLSE). 



This equation is given by 

iqt+qxx+'2q\q\'^+'ji {qxxxx + 6g^g* + 4g|g^p + 8|gpg^^ + 2g^g*^ + 6|g|^g) = 0. 

(1) 
Here q{x,t) is the complex envelope and 71 denotes the strength of higher 

order linear and nonlinear effects. When we consider the propagation of ultra 

short pulse propagation through optical fibre, i.e, less than 100 fs pulses, it 

has been shown that higher order dispersion, self steepening, self frequency 

and quintic effects should be included in the model. The above equation 

has been shown to be integrable and admits exact soliton solutions and also 

gauge equivalent to Heisenberg spin chain equation. Thus it is an interesting 

problem to find how these higher order terms will affect the breather and 

rogue wave in an associated optical system/ spin system by means of changing 

the value of 71. It is our prime aim to answer this problem in this paper. 

The paper is organized as follows. In Section 2, the Lax pair and the 

Darboux transformation (DT) are introduced. In Section 3, we shall give 

the first order breather and it's limit of infinitely large period. In Section 4, 

higher order rogue waves are given. The compression effects on the breather 

and rogue waves produced by the higher order terms of GNLSE are given in 

sections 3 and 4 respectively. Section 5 is devoted to conclusions. 

2. Lax Pair and Darboux transformation 

As discussed above, in this section, we would like to recall the Lax pair 
of GNLSE [SH ESI E3|and to show its Darboux transformation. According to 
the AKNS formalism. Lax pair for eq.rtT| is written as 

$^ = M$,<l>i = iV$ (2) 

with the following matrixes 

M = tXUo + f/i = zA ( "q^ ^i)^{ -q* 
and A^ = 8^71^4 - 2iV2. Here 

y _ f A2 - 1 qq* iqX - \q^ \ _ 
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M =A^ - 2^^*^^ + 4 ^^^*x - (1^(1*) ^ + g (3?^?*^ + q* qx,x + qql,^ - Qxql), 



Bi =iqX - -q^X 



Ci = - iq* X^ 
Moreover, 



^ T {qx,x + '^q^q*) X + -{qx,x,x + 6 qq* qx) 
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^q*x>^^ + 4 {q*x,x + 2gg*^) A + -(g*,^,^ + 6gg* ql). 
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0(A) \ ^ /0(A;x,t) 
^(A) ; VV^(A;x,t) 

denotes the eigenfunction of Lax pair eq.([2]) associated with A. 

The Lax pair of GNLSE provides a basis for the solvabihty of this equa- 
tion by means of the Darboux transformation. To construct the n-fold DT, it 

is necessary to introduce 2n eigenfunctions /«=(/' ) = $(Aj) associated 

with eigenfunction Xi{i = 1,2, ■■■ ,2n) and satisfy corresponding reduction 
condition X2k = ^2k-i ^^ ^^ have done for the NLS equation [311 ES]- Further- 
more, a similar n-fold DT determinant representation derived for the NLS 
equation as given in [SH? ] For example, from the one-fold DT, we get 



rW 



with Ai = ^+ir], X2 = XI = ^-if], /2,i 



2iAi 

"/r,2)/2,2 = /l*l7 ^1 



(3) 
/i,i Ai/1^1 

/2,1 A2/2,l 
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/l,l /l,2 
/2,1 /2,2 



/l,l/2,2 — /l, 2/2,1- 



3. The first order breather and its limit 

In this section, we first solve the eigenfunctions associated with a periodic 
seed q^^\ and then use it to get a first order breather by using the determinant 
representation of one-fold DT in eq.([3]). Further, this breather implies a first 
order rogue wave in the limit of infinitely large period. These two solutions 
have 71 explicitly such that we can use it to study the effects of breather and 
rogue waves affected by the higher order terms. 

Considering a periodic solution in following form 



ce 



«p 



(4) 



with p = ax + bt,b = (a^ — 12a^c^ + 6c^) 71 + 2 c^ — a^. By the method of 
separation of variables and the superposition principle, we have the following 
eigenfunction associated with g'"' 

h,,= hcei'^^')+^h[^ + h + x)ei'^-'), (5) 

fi,2= hce-<'^^') + zhl^ + h + X) e-i'^-') , (6) 



with, h=^c^+{\ + ^Y = hn + ihiM = e""^"' "^"' ''\k2= e-'^'i'' '+'' ^') , 

d =-fi (a(a^-6c^) -8A^ + 4aA^+ {Ac^ - 2 a^) \) ht + {x + {2\ - a) t) h 
= (x+ (2A-a + 7i {a{a^ -6c^) -8A^ + 4aA^+ {A c^ - 2 a^) \)) t) h 
= {x + {dji + idi)t)h 

Using the one-fold DT,a first order breather is constructed in the form 

[^] _/ 27/ { [wi cos (2G) - W2 cosh (2F)] -i[{wi- 2c^) sin (2G) - W3 sinh (2F)] } \ ^^ 
^ "^"^"^ wi cosh (2F) - W2 cos (2G) j '^ ' 

(7) 

with wi = c^ + {hi + rif + {i + hR+ f)2,W2 = 2c(/i/ + vi),w^ = 2c(^ + 

hR + |), -F = x/i/ + {dnhj + dihfi)t,G = xhji + {dj^hji — djhj)t. This is 

a periodic traveling wave. The coefficient 71 can affect the period of the 

breather through G. 

It is trivial to find 

|gWp(0,0) = (c + 2r/f, 

which is the height of peaks of this breather. Obviously, the height is in- 
dependent of a,^ and 71. This does not mean that 71 cannot affect the 
properties of the breather. In fact, we can see from eq.([7]) that 71 actually 
controls the period of the breather. This observation can be clearly seen in 
Figure [T] the number of peaks on same time interval is increasing when 71 
goes up from to 3 with a constant gap 0.5. These pictures clearly show 
that the resulting breather is compressed by the higher order effects due to 
the presence of 71 than 71 = case. In addition to the above, when the value 
of 7i increases, the number of peaks also increases. We use a short interval 
in Figure fIFa) and (b) to avoid too many peaks in it. 

Now we can consider what will happen in a breather when its period goes 
to infinity. According to the explicit expression in eq.([7]) of the first order 



breather, this hmit can be reahzed by setting X = ^ + ir] -^ \q = —- + ic, 

i.e., hm^^Ao Q^^^- For siniphcity, on setting a = , the hmit of the breather 
solution is obtained as 

with T = 4 (1 + 6 7iC^) c^t, X = 2cx. This is nothing but a first order rogue 
wave possessing asymptotic height 1 when x and t go to infinity. Further we 
find iQiimitlmaxi^^^) = 9^^' which denotes the height of a first order RW. 
We can see from T that,hke in the case of breather compression, 71 also 
responsible for compression effect of RW in the time direction, which is clearly 
seen in Figure [2] with 7i = 0,0.5,l,1.5 respectively. 

4. Higher order Rogue Waves 



The limit method in eq. rtSj) is not applicable for the higher order breather 
when Aj — > Ao(^ > 2). We can overcome this problem by using the coefficient 
of the Taylor expansion in the determinant representation of a higher order 
breather gt"! [5^ |55]. Similar to the case of NLS |50j, the first order rogue 
wave of GNLSE is given by 



Ql^l = [^^J^ - Ij ^e-. (9) 

Here Fi =4, Gi = 16 (1 + 6 71 c^ - 6 0^71) cH, 

H^ =4 c^x^ + (32 c^-iia^ - 16 c^ (l + 12 71 c^) a) xt + (64 c^ji'^a'^ 
- 64 0^71 (1 + 3 7ic2) ^4 ^ ^g ^2 (^^ ^ ^2 7ic2 ^ 72 cV^) a^ 

+ 16c^(67ic2 + l)V + l. 

Note that gl-J, reduces to the gjjj^jj when a = 0. 

What follows is the second-order rogue wave given by the Taylor expan- 
sion when Aj — )■ Ao(« = 1, 3) as the case of NLSfnO]- There are two patterns 
for the second order RW. The first one is called the fundamental pattern 
possessing a highest peak surrounded by four small equal peaks in two sides. 
Setting a = 0, Si = 0, then the Taylor expansion in the determinant of g^^^ 
provides 

q?l = { ^^^^"h!^"^ - 1) ce^^^^(^--^-^^)*. (10) 



Here 



F2 =5T^ +{QX^ + 34) T^ - QATtc^ + X^ + QX'^ - 3, 

G2 =T^ + (2 X^ + 10) T^ - 32 tc^T^ + (X^ - 14 X^ - 23) T + 32 tc^ (X^ + l) 

H2 =T^ + (3 X^ + 43) T^ - 64 tc^T^ + (3 X^ - 66 X^ + 307) T^ 

+ (192 c^tX^ - 1088 tc^) T + 1024 cH^ + X^ + 3 X^ + 27 X^ + 9, 



and X, T are defined in eq.rtSj). By comparing two cases with 71 = and 
7i = 1 in Figure [3} the compression effect in t direction is shown clearly. The 
second is a triangular pattern, which consists of three equal peaks. Setting 
S2 = 0,a = 0, si = 50 — 502, then an explicit formula of this pattern is 

^[2] _ f. 12(F2trig + iG2trig) \ i2cH3-fic^+l)t /-.-.N 

Hrwtrig ~ \ ^ tt I Ce ^ ' , yii.) 

\ J^2trig J 

with 

F2trig = 5 T^ + (24 (1 + 6 7ic2) ^^2^2 ^ 24 (34 71 c^ + 3) tc^ + 1200 c^) T 
-3 + 1200c^X + QX^ + X\ 

G2tri9 = T^ + (8 (1 + 6 71 c') tc'X^ + 8 (1 + 30 71 c") tc^ + 600 c") T" 
+4 (1 + 67ic^) tc^X^ + (-24 (1 + 147ic2) cH - 600 c^) X^ 
+4800 (1 + 6 71 c^) c^tX - 12 (5 + 46 71 c^) ch - 600 c^ 

^2tH9 = T^ + 3 X^T^ + (12 (86 71 c^ + 9) tc^ + 1200 c^) T^ 

+ (3 X^ + 3600 c^X) T^ - (3600 c^ + 72 (22 71 c^ + l) tc^) X^T 
+X^ + 3X^ + 27X2 + 14400 c^ (34710^ + 3) t 
+144 (11 + 228710^ + 1228 71^0^) fc^ + 720000 c^ 
-1200 X^c^ + 3600 c^X + 9, 

and X,T are defined in eq.(8). Figure 4 is plotted for the Iql^i^trigl'^ ^^ show its 
compression effect. Because of the explicit appearance oft in F2, G2, -f^2, F2trigi 

[21 [21 

G2trig, H2trig, Qrw and ql-wtrig ^^^ cau uot be able to derive from the corre- 
sponding RWs of the NLSE by a mere scalar transformation of x, t. 

Next, we construct the third order RW of GNLSE by substituting present 

/j(z = 1, 2, ■ ■ ■ ,6) in gt3][^5Qj_ There exists are three patterns: a fundamental 

[31 [31 

pattern g™ with Si = S2 = 0, a triangular pattern ql^^rig with a = S2 = 

and a circular pattern ql^^irc with a = Si = 0. The explicit form qlw of 



the fundamental pattern with the third order RW is given in appendix and, 
for brevity, the very lengthy forms of other two cases are deleted. Note 

[31 

that qrw includes the t dependence explicitly. This fact shows as the second 
order RW solution cannot be obtained from the corresponding solution of 
the NLS equation by a scalar transformation of x, t. Furthermore, to show 
the compression effect on the RWs, Figures (|5|J7]) are plotted for different 
parametric choices. 

5. Conclusions 

In this paper, we considered the integrable version of the generalized non- 
linear Schrodinger equation with several higher order nonlinear terms which 
describes ultra short pulse propagation through nonlinear silica fibre and 
soliton type nonlinear excitations in classical Heisenberg spin chain. Using 
Daurboux transformation and periodic seed solutions, we have constructed 
the first order breather solution and also discussed the behavior of these so- 
lutions with an infinitely large period. Finally, we have also constructed the 
first order, second order and third order rogue wave solutions by the Taylor 
expansion. All of these solutions have parameter 71 denoting the contribution 
of higher order nonlinear terms. The compressed effects of these solutions 
are discussed through numerical plots by increasing the value of 71. This 
new phenomenon of the rogue wave is useful for us to observe or analyze 
its evolution in some complicated physical system. Another advantage of 
our results of this paper is that, as GNLSE is equivalent to spin chain, the 
rogue wave nature of spin systems can also be explained through suitable 
geometrical/ gauge equivalence method. In addition, as higher order linear 
and nonlinear effects in optical fibres are playing key role in explaining the 
generation and propagation of ultra short pulse through a silica wave guides, 
we hope that our results with all these higher order effects can be observed 
in real experiments in the near future. 

Acknowledgemnts This work is supported by the NSF of China under 
Grant No.l0971109 and 11271210, K.C.Wong Magna Fund in Ningbo Uni- 
versity. Jingsong He is also supported by the Natural Science Foundation of 
Ningbo under Grant No.2011A610179. Jingsong He thanks Prof. Fokas for 
his support during his visit at Cambridge from January to December 2012. 
KP wishes to thank the DST, DAE-BRN and CSIR, Government of India, for 
the financial support through major projects. Lihong Wang is also supported 
by the Natural Science Foundation of China under Grant No. 11074136 and 

9 



11101230, the Natural Science Foundation of Zhejiang province under Grant 
No. 2011R09025-06. 

References 

[1] G. P. Agrawal, Nonlinear Fibre Optics, Academic Press, San Diego 
(2006). 

[2] A. Hasegawa, M. Matsumoto, Optical solitons in fibres, Springer- 
Verlag, Berlin, (2003). 

[3] L.F. Mollenauer, J. P. Gordon, Solitons in Optical Fibres: Fundamen- 
tals And Applications, Academic Press, New York, (2006). 

[4] K. Porsezian and V. C. Kuriakose, Optical Solitons: Theory and Ex- 
periment, Lecture Notes in Physics, Vol.613, Springer- Verlag, Berlin, 
February (2003). 

[5] K. Porsezian and V. G. Kuriakose, Solitons in Nonlinear Optics: Ad- 
vances and Applications, European Journal of Physics-Special Topics, 
Springer- Verlag, Berlin (2009). 

[6] J. R. Taylor, Optical Solitons: Theory and Experiment, (Gambridge 
University Press, 1992). 

[7] F. Abdullaev, S.Darmanyan, P.Khabibullaev, Optical Solitons, 
(Springer Series in Nonlinear Dynamics, Berlin 1993). 

[8] F. Abdullaev, Theory of Solitons in Inhomogeneous Media, (John Wi- 
ley, New York 1994). 

[9] Y. S. Kivshar, G.P. Agrawal, Optical Solitons: From Fibres to Photonic 
Grystals, (San Diego: Academic, 2003). 

[10] G. P. Agrawal, Applications of Nonlinear Fibre Optics, Second Edition, 
Academic Press, New York (2006). 

[11] N. N. Akhmediev and A. Ankiewicz, Solitons, Ghapman Hall, London 
(1997). 

[12] D.H. Peregrine, J. Austral. Math. Soc. B 25(1983), 16-43. 

10 



[13] N. Akhmediev, A. Ankiewicz, M. Taki,Phys. Lett. A 373(2009), 675- 
678. 

[14] D. R. Solli, C. Ropers, P. Koonath,B. Jalali, Nature (London) 
450(2007), 1054-1058. 

[15] D. R. Solli, C. Ropers, B. Jalali,Phys. Rev. Lett.l01(2008),233902. 

[16] J. M. Dudley, G. Genty and S. Coen, Rev. Mod. Phys. 78(2006), 1135; J. 
M. Dudley, G. Genty, F. Dias, B. Kibler, N. Akhmediev, Opt. Express 
17(2009), 21497. 

[17] A. N. Ganshin, V. B. Efimov, G.V. Kolmakov, L. P. Mezhov-Deglin, 
P.V. E. McClintock,Phys.Rev.Lett.l01(2008),065303. 

[18] Yu. V. Bludov, V. V. Konotop,N. Akhmediev,Phys.Rev.A.80 
(2009),033610. 

[19] M.S. Ruderman,Eur. Phys. J. Special Topics 185(2010), 57-66. 

[20] W.M.Moslem, P.K.Shukla, B.Eliasson, Euro. Phys. Lett. 96 
(2011),25002. 

[21] R. Hohmann, U. Kuh, H.-J. Stockmann, L. Kaplan, E. J. Heller, 
Phys.Rev.Lett.l04(2010),093901. 

[22] M. Shats, H. Punzmann, H. Xia, Phys.Rev.Lett. 104(2010), 104503. 

[23] S.Vergeles,S.K.Turitsyn,Phys.Rev.E 83(2011), 061801 (R). 

[24] F. T. Arecchi, U. Bortolozzo, A. Montina, S. Residori, Phys.Rev.Lett. 
106(2011), 153901 

[25] A.Chabchoub, N.P.Hoffmann, N.Akhmediev, Phys. Rev. Lett. 106 
(2011),204502. 

[26] A. Chabchoub, N. P. Hoffmann, N. Akhmediev, J. Geophys. Res. 
117(2012), C00J02. 

[27] B. Kibler, J. Fatome, C. Finot, G. Millot, F. Dias, G. Genty,N. Akhme- 
diev, J. M. Dudley, Nature Phys.6(2010), 790-795. 



11 



[28] A. Ankiewicz, J. M. Soto-Crespo, N. Akhinediev, Phys. Rev. E. 
81(2010), 046602. 

[29] Y.S.Tao ,J.S. He,Phys. Rev. E. 85(2012),026601. 

[30] G.G.Yang, L.Li, S.T.Jia, Phys.Rev.E 85(2012),046608. 

[31] S.W.Xu, J.S. He, L.H.Wang, J. Phys. A 44(2011), 305203 

[32] S.W.Xu, J.S.He,J. Math. Phys.53(2012), 063507. 

[33] J.S.He, S.W.Xu, K. Porsezian, J. Phys. Soc. Jpn. 81(2012), 124007. 

[34] J.S. He, S.W. Xu, K. Porsezian, J. Phys. Soc. Jpn. 81(2012),033002. 

[35] C.Z.Li,J.S.He, K.Porsezian, Phys. Rev. E87(2013),012913. 

[36] U. Bandelow, N. Akhmediev, Phys. Rev. E86 (2012), 026606. 

[37] A. Ankiewicz, N. Akhmediev, J. M. Soto-Crespo, Phys. Rev. 
E82(2010),026602. 

[38] Yu. V.Bludov, V.V.Konotop, N. Akhmediev, Eur. J. Phys. Special Topics 
185(2010), 169-180. 

[39] B.L.Guo, L.M.Lin, Chin.Phys.Lett. 28(2011), 110202. 

[40] F. Baronio, A. Degasperis, M. Confortil,S. Wabnitz, Phys. Rev. Lett. 
109(2012),044102. 

[41] Z.Y. Qin, M.Gu, Phys.Rev.E 86(2012),036601. 

[42] A. Ankiewicz, N. Devine, N. Akhmediev, Phys.Lett.A373(2009),3997- 
4000. 

[43] M.Taki,A.Mussot, A. Kudhnski, E. Louvergneaux, M. Kolobov, M. 
Douay, Phys. Lett. A 374(2010), 691-695. 

[44] Z.Y. Yan, Phys. Lett. A. 374(2010), 672-679. 

[45] Y.Y.Wang, J.S.He,Y.S. Li, Commun.Theor. Phys. 56(2011), 995-104. 

[46] L.Wen, L.Li, Z.D.Li, S.W.Song, X.F.Zhang,W.M.Liu, Eur. Phys. J. D 

64 (2011), 473-478. 

12 



[47] S. W. Xu, J. S. He, L. H. Wang, Europhys. Lett.97(2012), 30007. 

[48] C. Q. Dai, G.Q. Zhou,J.F. Zhang, Phys. Rev. E 85(2012), 016603. 

[49] J. M. Soto-Crespo, Ph. Grelu, N. Akhmediev, Phys.Rev. E84(2011), 
016604. 

[50] J.S.He, H.R. Zhang, L.H. Wang, K. Porsezian, A.S.Fokas,A generating 
mechanism for higher order rogue waves ( jar Xiv: 1209.3742^ 3) . 

[51] M.Lakshmanan, K.Porsezian, M.Daniel, Phys. Lett. A,133(1988), 483- 
488. 

[52] K. Porsezian, M. Daniel, M. Lakshmanan, J. Math. Phys. 33(1992), 
1807-1816. 

[53] K. Porsezian, Phys. Rev. E 55(1997), 3785-3789. 

[54] V. B. Matveev, M. A. Salle, Darboux transformations and soli- 
tons (Springer, Berlin, 1991). 

[55] J. S. He, L. Zhang, Y. Cheng, Y. S. Li, Sci. China A 12(2006), 1867- 
1878. 



Appendix 



,E^( '^^\+^^^) +l)ce^-^(^^--^+^> 



F3 ^IIT^" + A5X^T^ + 180 {130 -/ic^ + 11) tc^T'^ + {70 X"^ + 420 X^) T^ 

+ (50 X^ + 480 (15300 71^0* + 2508 71 c^ + 73) t'^c^) T^ - 600 (38 71 c^ + l) tc^X^T^ 

+ (15 X^ + 7200 (2460 -fi^c'^ + 308 71 c^ + 15) t^c'^X^) T^ 

+ (-240 (38 7ic2 + 1) tc^X^ - 28800 (4056 €^71^ + 17 + 300-fi'^c^ + 210 71 c^) t^c^) T 

+X^° + 15X^ + 210X'^ + (-7200 (2207120"* + 20 71 c^- l) t^c** - 450) X'' 

+ (-43200 (628 7120-* + 124 71 c^ + 5) t^c^ - 675) X"^ + 675 

+10800 (2452 7120* + 2871 c^ - 3) t^c^ 
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Ga =T" + SX^T^ + 20 (102 71 c^ + 5) tc^T^ + {lOX^ - 60 X^) T^ 

+ (10X^+480 {l2-ii^c^ -2Q8-fic^ -29)fc^)T^ - 120 (82 71 c^ + 7) tc^X^T-* 

+ (5X* + 1440 (4524 7i2c'' + 548 71 c^ + 19) t^c^X"^) T^ + ( - 80 (907ic2 + 7) tc^X^ 

-5760 (99432 c<^7i3 + 28676 7i2c4 + 3086 71 c^ + 107) t^c^)T'^ 

+ (^-7200 (I47ic2 + iffc^X'^ + X^^\ T - 60 (14710^ + l) tc^X^ 

+120 (2 7ic2 - 5) tc^X^ - 1800 (IO71C2 + 3) tc^X^ 

+ (57600 (-126 71 c^ + 1176 0671^ - 5647i2c4 - 7) t^c^ + 2700 (170710^ + 7) tc^) X"^ 

+18900 (I47ic2 + 1) tc^ - 14400 (ll + 1254 71 c^ + 1808471^0'' + 84168 c^7i^) t^c^ 

H3 =T^^ + 6X^T^° + 24 (206 71 c^ + 2l) tc^T^ + (15 X'^ + 270 X^) T^ 

+ (20 X*^ + 720 (7596 7i2c'' + 1636 71 c^ + 83) t^c^) T^ - 240 (42 7ic2 - l) tc^X^T^ 

+ (15 X^ + 8640 (3012 ji^c^ + 492 71 c^ + 25) t^c^X'^) T'^ 

+ (-240 (82 7ic2 + 3) tc^X^ - 57600 (3048 0^71^ - 2604 7i2c'' - 450 71 c^ - 17) t^c^) T^ 

+ (6^1° - 21600 (5487i2c'' + 76710^ + l) t^c^X^) T^ 

+ (172800 (46968 c^7i3 + 11196 7i2c'* + 906 71 c^ + 29) t^c^X"^ 

-360 {22"/ic^ + l)tc^X^)T 

+X^'^ + 6X^° + 135 X^ + (2880 (13247120"* + 164 71 c^ + 3) fc^ + 2340) X^ 

+ (-43200 (42871^0'' - 12 7ic2 - 5) t^c'^ + 3375) X^ 

+ (-64800 (25007i2c'* +492 7ic2 +9) t^c^ + 12150) X^ + 2025 

+64800 (7260 71^0'' + 772 71 c^ + 23) t^c'^ 

+172800 (213 + 8056 71 c^ + 720096 0^71^ + 115128 7i2c'* + 1836624 71 ''c^) t^c^ 
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(a) a = ^ = 0,c= 1,77= i 71 = 1 



(b) a = e = 0,c=|,7/= 1,71 = 1 




(c) a = ^ = 0,c= 1,77= i,7i = 1 




(d) a = C = 0,c= |,7?= i 71 = 





(e)a = e = 0,c=|,77=i,7i = 2 (f) a = ^ = 0,c = §,77 = ^,71 = 3 

Figure 1: (Color online) The dynamical evolution of the first order breather jg'^'P on the 
(a;, i)-plane. When the value of 71 increases, the number of peaks on same interval of time 
also increases. 
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(a) c= i,7i = 




(b) c= i,7i == i 




(c) c= i,7i = 1 




(d)c=i,7i = 2 



Figure 2: (Color online) The dynamical evolution of the first order rogue wave Igjjj^jtP 
on the (x, i)-plane. For larger values of 71, it is clear that the compression effects in t 
direction are quite high. 
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(a) a = si =0,c= ^,7i =0 




(b)a: 



0,c=^,7i=0 




(c) a = si = 0, c= ^,7i = 1 



i^r 6- 




^/2' 



(d) a = Si = 0, c = ^, 7i = 1 



%/2' 



Figure 3: (Color online) The dynamical evolution of the second order rogue wave |qru;p 
on the (x, t)-plane. Comparing of (a) and (b) with (c) and (d) indicates effective high 
compression in the t direction. 
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-6 -I -2 a 




(a) a == 0, c = ;35, si = 50- (b) o = 0, c = 3^' ^i = ^0 - 50i, 71 = 
50i,7i =0 ' 



V2 



72' 





(c) a = 0, c= 45,si = 50-50i,7i = ^ (d) a = 0, c = ^, si = 50-50i,7i = 



v^ 



\/2' 



Figure 4: (Color online) The dynamical evolution of the second order rogue wave l^riitrioP 
on the {x, t)-plane. It is shown from (b),(c) and (d) that rogue wave compression increases 
as the value of 71 increases. 




(a) a == si = 0,c= ^,71 = (b) a = si = 0, c = ^,71 = 



72 



72' 




(c) a = si =0,c= ;^,7i = i 




(d) a = si = 0, c = ^ 



72 



71 



Figure 5: (color online) The dynamical evolution of the third order rogue wave Iq™^ on 
the (a;, i)-plane. By comparison with (a) and (b), (c) and (d) are highly compressed. 
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(a) c 




72 



, Si — — 50i,7i — 



(b) c=^,si=-50i,7i=0 




(c)c: 



v^ 



,Sl 



-50i,7i 




(d)c 



V2 



, Si 



-50i,7i 



Figure 6: (Color onlinc)Thc dynamical evolution of the third order rogue wave |9™triQp 
on the (a;,i)-plane. It is shown from (b),(c) and (d) that rogue wave is compressed more 
while increasing the value of 71 . 



20 





(a) a = 0,c= -^,S2 = 5000i,7i =0 (b) a = 0, c = ^, 53 = 5000i,7i =0 





(c) a = 0,c= ^,S2 = 5000i,7i = I (d) a = 0, c = ^, 53 - 5000i,7i = 1 



irwcircl 



Figure 7: (Color online) The dynamical evolution of the third order rogue wave \q^ 

on the (x,i)-plane. It is shown from (b),(c) and (d) that rogue wave is changing its shape 

and also compression increases by increasing the value of 71 . 
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